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Understanding the consequences of the gravitational interaction between a star and a planet 
is fundamental to the study of exoplanets. The solution of the two-body problem shows that 
the planet moves in an elliptical path around the star and that each body moves in an ellipse 
about the common center of mass. The basic properties of such a system are derived from first 
principles and described in the context of detecting exoplanets. 



1. INTRODUCTION 

The motion of a planet around a star can be understood in 
the context of the two-body problem, where two bodies ex- 
ert a mutual gravitational effect on each other. The solution 
to the problem was first presented by Isaac Newton (1687) 
in his Principia. He was able to show that the observed el- 
liptical path of a planet and the empirical laws of planetary 
motion derived by Kepler (1609, 1619) were a natural con- 
sequence of an inverse square law of force acting between 
a planet and the Sun. According to Newton's universal law 
of gravitation, the magnitude of the force between any two 
masses mi and m 2 separated by a distance r is given by 



F = G 



m\m,2 



(1) 



where G = 6.67260 x 10~ n Nm 2 kg~ 2 is the universal 
gravitational constant. The law is applicable in a wide va- 
riety of circumstances. For example, the two bodies could 
be a moon orbiting a planet or a planet orbiting a star. New- 
ton's achievement was to show that motion in an ellipse is 
the natural consequence of such a law. A more difficult 
task is to find the position and velocity of an object in the 
two-body problem; this is commonly referred to as Kepler 
problem. In this chapter we derive the basic equations of 
the two-body problem and solve them to show how ellip- 
tical motion arises. We then proceed to solve the Kepler 
problem showing how motion around the common center 
of mass of the two-body system can be used to infer the 
presence of planetary companions to a star. Finally we give 
a few representative examples among extra-solar planets al- 
ready detected. For the most part we follow the approach of 
Murray & Dermott (1999). 



2. BASIC EQUATIONS 

Consider a star and a planet of mass mi and m^, respec- 
tively, with position vectors ri and r 2 referred to an origin 
O fixed in inertial space (Fig. 1). 
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Fig. 1. — The forces acting on a star of mass mi and a planet of 
mass m,2 with position vectors ri and r-2. 



The relative motion of the planet with respect to the star 
is given by the vector r = T2 — r%. The gravitational forces 
acting on the star and the planet are 



Fi = mii 1 ! 



mim 2 

"(_t 5 r . 



F 2 = m 2 r 2 = -G 5 — r 



(2) 
(3) 



respectively. Now consider the motion of the planet m 2 
with respect to the star mi. If we write r = r 2 — ?i we can 
use Eq. (1) to obtain 



r + G(m\ + m 2 W = . 



(4) 



If we take the vector product of r with Eq. (4) we have 
r x f = which can be integrated directly to give 



r x r = h 



(5) 



where h is a constant vector which is simultaneously per- 
pendicular to both r and r. Therefore the motion of the 
planet about the star lies in a plane (the orbit plane) perpen- 
dicular to the direction defined by h. Another consequence 
of this result is that the position and velocity vectors will 
always lie in the same plane (see Fig. 2). Equation (5) is 
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Fig. 2. — The motion of m,2 with respect to mi defines an orbital 
plane (shaded region), because r x r is a constant vector, h, the 
angular momentum vector and this is always perpendicular to the 
orbit plane. 
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Fig. 3. — The geometry of the ellipse of semi-major axis a, semi- 
minor axis b, eccentricity e and longitude of periapse zu. 



often referred to as the angular momentum integral and h 
represents a constant of the two-body motion. 

In order to solve Eq. (4) we transform to a polar coor- 
dinate system (r, 9) referred to an origin centered on the 
star with an arbitrary reference line corresponding to 9 = 0. 
In polar coordinates the position, velocity and acceleration 
vectors can be written as 



;• r 

r f 



r6\ 



(r - r9 2 )r + 



rdt V 



(6) 
(7) 

(8) 



where f and 9 denote unit vectors along and perpendicular 
to the radius vector respectively. Substituting Eq. (7) into 
Eq. (5) gives h = r 2 9 z, where z is a unit vector perpendic- 
ular to the plane of the orbit forming a right-handed triad 
with f and 9. The magnitude of this vector gives us 

h = r 2 9. (9) 

Therefore, although r and 9 vary as the planet moves around 
the star, the quantity r 2 9 remains constant. The area ele- 
ment dA swept out by the star-planet radius vector in the 
time interval dt is given in polar coordinate by 



dA 



f 

Jo 



r dr d9 



1 



l d9. 



and thus 



A = —r 2 9 = —h = constant. 
2 2 



(10) 



(11) 



This is equivalent to Kepler's second law of planetary mo- 
tion which states that the star-planet line sweeps out equal 
areas in equal times. 

Using Eq. (6) and comparing the f components of 
Eqs. (4) and (8) gives the scalar differential equation 



r9 2 



G(mi 



m 2 ) 



(12) 



In order to find r as a function of 9 we need to make the 
substitution u = 1/r. By differentiating r with respect to 
time and making use of Eq. (9) we can eliminate time in the 
differential equation. We obtain 



h d2u a h2 2 d2 « 
r = ~h W e = -hu — 

and hence Eq. (12) can be written 

d 2 u G(m,i + rn 2 ) 

d6 2+U = h? • 



(13) 



(14) 



This is a second order, linear differential equation, often 
referred as Binet's equation, with a general solution 



G(mi + m 2 ) , 
/i 2 [ 



w)} , 



(15) 



where e (an amplitude) and m (a phase) are two constants 
of integration. Substituting back for r gives 



P 



1 



ecos 



(16) 



where p = h 2 /G(rrii + 1712). This is the general equation 
in polar coordinates of a set of curves known as conic sec- 
tions where e is the eccentricity and p is a constant called 
the semilatus rectum. For a given system the initial con- 
ditions will determine the particular conic section (circle, 
ellipse, parabola, or hyperbola) the planet follows. We con- 
sider only elliptical motion for which 



p = a(l - e 2 ) , 



(17) 



where a, a constant, is the semi-major axis of the ellipse. 
The quantities a and e are related by 



b 2 = a 2 (\-e 2 ), 



(18) 



where b is the semi-minor axis of the ellipse (see Fig. 3). 

Therefore, for any given value of 9 the radius is calcu- 
lated using the equation 



a(l-e 2 ) 



(19) 
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Hence the path of the planet around the star is an ellipse 
with the star at one focus; this is Kepler's first law of plan- 
etary motion. Note that in the special case where e = (a 
circular orbit), r = a and the angle w is undefined. 

The angle 8 is called the true longitude. Equation (19) 
shows that the minimum and maximum values of r are 
a(l — e) (at = w) and a(l + e) (at 6 = w + n), re- 
spectively. These points are referred to as the periapse and 
the apoapse, respectively, although for motion around a star 
they can also be referred to as the periastron and apastron. 

The angle w (pronounced "curly pi") is called the lon- 
gitude of periapse or longitude of periastron of the planet's 
orbit and gives the angular location of the closest approach 
with respect to the reference direction. If we define the true 
anomaly to be the angle / = 6 — ru (see Fig. 3) then / is 
measured with respect to the periapse direction and Eq. (19) 
can be written 



a(l 



1 + e cos / 



(20) 



In this case if we define a cartesian coordinate system cen- 
tered on the star with the x-axis pointing towards the peri- 
apse (see Fig. 3), then the position vector has components 



x = r cos / 
y = r sin / . 



(21) 
(22) 



Although we have eliminated the time from our equation 
of motion, we can relate the orbital period, T, to the semi- 
major axis, a. The area of an ellipse is A — nab and this is 
swept out by the star-planet line in a time, T. Hence, from 
Eq. (11), A = hT/2nnd so 



4tt 2 



G(mi + m 2 ) 



(23) 



This is Kepler's third law of planetary motion. It implies 
that the period of the planet's orbit is independent of e and 
is purely a function of the sum of the masses and a. If we 
define the mean motion, n of the planet's motion as 



2tt 



then we can write 



2„3 



G(m\ + m 2 ) = n a 



(24) 



(25) 



and hence 



h 



na 



V 1 - e 2 = v/GOi + m 2 )a(l - e 2 ) . (26) 



There is an additional constant of the two-body motion 
which is useful in calculating the velocity of the planet. 
Taking the scalar product of f with Eq. (4) and using 
Eqs. (6) and (7) gives the scalar equation 



r • r 



+ G(mi + m 2 )- 







(27) 



which can be integrated to give 
1 



2 _ Gjmi +m 2 ) _ 

2v o , 

r 



(28) 



where v 2 = f • r is the square of the velocity and C is a 
constant of the motion. Equation (28) is called the vis viva 
integral. It shows that the orbital energy per unit mass of 
the system is conserved. 

Because w is a constant, 9 = f and Eq. (7) gives 



v = r • r = r + r j . 

By differentiating Eq. (20) we obtain 

r f e sin f 
r = — — 

1 + e cos / 

and hence, using Eqs. (9) and (16), we have 

na 



■ e sin / 



and 



rf 



na 



:(l + ecos/) . 



Therefore we can write Eq. (29) as 

2 2 

v 2 = ^-^2(l + 2ecos/ + e 2 ). 



(29) 



(30) 



(31) 



(32) 



(33) 



This shows the dependence of v on /. A little further ma- 
nipulation gives 

v 2 = G{m 1 +m 2 )(---) (34) 



which shows the dependence of v on r. 



3. SOLUTION OF THE KEPLER PROBLEM 

In the previous section we solved the equation of motion 
of the two-body problem to show the path of the planet's 
orbit with respect to the star. However, in the process we 
eliminated the time and so although we can calculate r for 
a given value of 9, we have no means of finding rasa func- 
tion of time. This is the essence of the Kepler problem. 

Our starting point is to derive an expression for f in terms 
of r. We can do this by using Eqs. (20), (32) and (34) to 
rewrite Eq. (29) as 



f 2 = n 2 a 3 {--- 



n 2 a 4 (l - e 2 ) 



This simplifies to give 



f = — \J a 2 e 2 — (r — a) 2 . 
r 



(35) 



(36) 
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In order to solve this differential equation we introduce 
a new variable, E, the eccentric anomaly, by means of the 
substitution 



r = a(l — ecosE) . 
The differential equation transforms to 

E = n - . 

1 — e cos E 

The solution can be written as 

n(t — t Q ) = E~esmE, 



(37) 



(38) 



(39) 



where we have taken t to be the constant of integration and 
used the boundary condition E = when t = to . At this 
point we can define a new quantity, M, the mean anomaly 
such that 



M = n(t - t ) • 



(40) 



where to is a constant called the time ofperiastron passage. 
There is no simple geometrical interpretation of M but we 
note that it has the dimensions of an angle and that it in- 
creases linearly with time. Furthermore, M = f = when 
t = t a or t = to + T (periapse passage) and M = f = it 
when t = t + T/2 (apoapse passage). We can write 



M = E - e sin E . 



(41) 



This is Kepler's equation and its solution is fundamental to 
the problem of finding the orbital position at a given time. 
For a particular time t we can (i) find M from Eq. (40), (ii) 
find E by solving Kepler's equation, Eq. (41), (iii) find r 
using Eq. (37), and finally / using Eq. (20). 

The key step is solving Kepler's equation and this is usu- 
ally done numerically. Danby (1988) gives several numeri- 
cal methods for its solution. For example, if we define the 
function 



g(E) = E -esmE - M 



(42) 



then we can use a Newton-Raphson method to find the root 
of the non-linear equation, g(E) = 0. The iteration scheme 
is 



Ei 



+i 



Ei 



9_m 



0,1,2,... 



(43) 



where g'(Ei) = dg(Ei)/dEi = 1 - e cos Ei and the iter- 
ations proceed until convergence is achieved. A reasonable 
initial value is Eq = M since E and M differ by a quantity 
of order e (see Eq. (41)). 

Although we cannot have an explicit relation between 
the angles / and M, from Kepler's second law (Eq. 9) and 
(Eq. 26) it is possible to write: 

&f = ny/l-e 2 (^) dt = y/ \ - e 2 AM . (44) 



The above relation is useful when we want to average any 
physical quantity over a complete orbit. For instance, 



i \ _ i r 2n dM _ i 

^2/ 2W ~ ~ a 2yr^P ' 



(45) 



To complete the set of useful angles we define the mean 
longitude, A by 



A = M + w . 



(46) 



Therefore A, like M, is a linear function of time. It is im- 
portant to note that all longitudes (8, w, A) are defined with 
respect to a common, arbitrary reference direction. 



4. THE ORBIT IN THREE DIMENSIONS 

Of the orbital elements we have defined so far, two (a 
and e) are related to the physical dimensions of the orbit 
and the remaining two (w and /) are related to orientation 
of the orbit or the location of the planet in its orbit. Note 
that there are many alternatives to / (e.g. 6, M and A) and 
the time of periapse passage, to, can be used instead of / 
since the latter can always be calculated from the former. 
We have already noted that w is the angular location of the 
periapse direction measured from a reference point on the 
orbit. 

Consider the planet's position vector, 



r = (x, y, 0) = xx + y y + z 



(47) 



in a three-dimensional coordinate system where the x-axis 
lies along the major (long) axis of the ellipse in the direction 
of periapse, the y-axis is perpendicular to the a>axis and lies 
in the orbital plane, while the z-axis is mutually perpendic- 
ular to both the x- and j/-axes forming a right-handed triad. 
By definition orbital motion is confined to the x-y plane. 
Consider a standard coordinate system where the direction 
of the reference line in the reference plane forms the X- 
axis. The F-axis is in the reference plane at right-angles 
to the X-axis, while the Z-axis is perpendicular to both the 
X- and Y"-axes forming a right-handed triad. 

Let I denote the inclination, the angle between the orbit 
plane and the reference plane. The line formed by the in- 
tersection of the two planes is called the line of nodes. The 
ascending node is the point in both planes where the or- 
bit crosses the reference plane moving from below to above 
he plane. The longitude of ascending node, ft is the angle 
between the reference line and the radius vector to the as- 
cending node. The angle between this same radius vector 
and the periapse of the orbit is called the argument of pe- 
riapse, u). Note that the inclination is always in the range 
< I < 180°. An orbit is said to be prograde if I < 90° 
while if I > 90° the motion is said to be retrograde. We 
can also define 



w = O + to 



(48) 
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Fig. 4. — The relationship between the (x, y, z) and (X, Y, Z) 
coordinate systems and the angles cj, I and fi. 



where vo is the longitude of periapse introduced above but 
that now, in general, the angles £1 and u> lie in different 
planes so that zu forms a 'dog-leg' angle. 

The orientation angles I, and 10 are illustrated in 
Fig. 5. It is clear that coordinates in the (x, y, z) system 
can be expressed in terms of the (X, Y, Z) system by means 
of a series of three rotations: (i) a rotation about the z-axis 
through an angle lj so that the x-axis coincides with the 
line of nodes, (ii) a rotation about the x-axis through an an- 
gle I so that the two planes are coincident and finally (iii) a 
rotation about the z-axis through an angle il. We can rep- 
resent these transformations by two 3x3 rotation matrices, 
denoted by P x ((f)) (rotation about the x-axis) and P z (</>) 
(rotation about the z-axis), with elements 



1 

cos <f> — sin < 
sin 6 cos < 



(49) 



and 



Consequently 



COS 


— sin 


sin <j) 


cos<? 









(50) 



and 



P,(n)P x (J)P,(w) I y 



y) =p-V)p- 1 (/)p- 1 (Q) ( y 



(51) 



(52) 



where P" 1 ^) = ~P x {-4>) and Pj 1 ^) = p ^(-0) are the 
inverse of the matrices of V x (<f)) and P z (</>), respectively. 
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Fig. 5. — The position vectors of star and planet with respect to 
the origin, O, and with respect to the center of mass of the star- 
planet system, O' . 



If we now restrict ourselves to coordinates which lie in 
the orbital plane, we have x — r cos f,y — r sin f,z = 
and 

X = r (cos cos(o; + /) — sin fl sin(oj + /) cos /) (53) 
Y = r (sin f2 cos(w + /) + cos ft sin(w + /) cos /) (54) 
Z =rsin(w-|-/)sin7. (55) 



5. BARYCENTRIC MOTION 

In order to determine the observable effects of an orbit- 
ing planet on a star it helps if we consider the motion in 
the center of mass or barycentric system (see Fig. 5). The 
position vector of the center of mass of the system is 



R 



toiTi + m 2 r 2 



mi + 77i 2 
From Eqs. (2) and (3) we have 

mif ! + m 2 r 2 



R 



mi + m 2 



= 0, 



(56) 



(57) 



and by direct integration R = V = constant. These equa- 
tions imply that either (i) the center of mass is stationary 
(the case when V = 0), or (ii) it is moving with a constant 
velocity (the case when V ^ 0) in a straight line with re- 
spect to the origin O. Then, if we write Ri = ri R and 
R 2 = r 2 — R, we have 



miRi + rn 2 R 2 = . 



(58) 



This implies that (i) Ri is always in the opposite direction 
to R 2 , and hence that (ii) the center of mass is always on 
the line joining mi and m 2 . Therefore we can write 



i?, 



R 2 = r. 



(59) 



where r is the separation of mi and m 2 , and the distances 
of the star and planet from their common center of mass are 
related by mi-Ri = — m 2 i? 2 (Eq. 58). Hence 



Ri = 



m 2 



mi + m 2 



and i?o = — - 



mi 



mi + m 2 



(60) 
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Fig. 6. — (a) The motion of the planet m-2 with respect to the 
star mi in the two-body problem; the dashed curve denotes the 
elliptical path of the center of mass, O' . (b) The motion of the 
masses mi and m-2 with respect to the center of mass, O' , for the 
same system. For the purposes of illustration we used mz/mi = 
0.2 ande = 0.5. 



Therefore each object will orbit the center of mass of the 
system in an ellipse with the same eccentricity but the semi- 
major axes is reduced in scale by a factor (see Fig. 6) 



a i 



m 2 



mi + m 2 



a and a 2 



mi + m 2 



(61) 



The orbital periods of the two objects must each be equal 
to T and therefore the two mean motions must also equal 
be equal (rii = n 2 — ri), although the semi-major axes are 
not. Each mass then moves on its own elliptical orbit with 
respect to the common center of mass, and the periapses of 
their orbits differ by ir (see Fig. 6b). 

We are now in a position to revisit the expression for the 
radial velocity of the star, v r . Observers usually take the 
reference plane (X, Y) to be the plane of the sky perpen- 
dicular to the line of sight, the Z-axis oriented towards the 
observer (Fig. 7). Thus, the radial velocity of the star is 
simply given by the projection of the velocity vector on the 
line of sight. Since ri = R + Ri this gives 



v r = ri • Z = Vz 



m 2 



mi + m 2 



(62) 



where Vz = V • Z is the proper motion of the barycenter 
and Z can be obtained directly from Eq. (55): 

Z = r sin(o; + /) sin I + rf cos(w + /) sin I , (63) 

or, making use of Eqs. (31) and (32), 

na sin / 



Z 



(cos(o; + /) + ecosw) 



We can now write 



v r — Vz + K (cos(w + /) + e cos oj) 



where 



K 



m 2 na sin I 



mi + m 2 y/l - e' 2 



(64) 



(65) 



(66) 
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Fig. 7. — The relationship between the star's velocity around the 
center of mass, ri, and its radial component along the line of sight, 



6. APPLICATION TO EXTRA-SOLAR PLANETS 

More than 500 extra-solar planets are known to date 1 , 
and the number is continuously rising. Looking at this data, 
we can admire the wide variety of possible orbital parame- 
ters and physical properties: central stars of spectral types 
from F to M, minimum masses from 2 Earth-masses to more 
than 20 times the mass of Jupiter, orbital periods of one day 
to more than fifteen years (the same time as the length of 
the observations), and eccentricities ranging from perfect 
circular orbits to extreme values of more than 0.9. There 
are planets as close as 0.014 AU and as far as 670 AU from 
their host stars. 

In Table 1 we report two examples of extreme values 
for the eccentricity, obtained using the radial velocity tech- 
nique. The first example (HD 156846 b) corresponds to a 
highly eccentric orbit, while the second one (HD 83443 b) 
shows an almost circular orbit. At present, both planets are 
in single-planet systems, which allows us to apply directly 
the formulae derived in previous sections to the analysis of 
their motion. 

6.1. HD 156846 b 

HD 156846 has been observed with the CORALIE spec- 
trograph at La Silla Observatory (ESO) from May 2003 to 
September 2007. Altogether, 64 radial velocity measure- 
ments with a mean uncertainty of 2.8 m/s were gathered. 
Figure 8a shows the CORALIE radial velocities and the cor- 
responding best-fit Keplerian model. The resulting orbital 



1 The Extrasolar Planets Encyclopedia, http://exoplanet.eu/ 
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Table 1 

two examples of extreme orbital eccentricities for extra-solar planets. 





HD 156846 b 


HD 83443 b 


discovery year 


zUU/ 


ZUUZ 


Utlltl I CI. 


Tamii7 <3f /// r?008^ 

1 ell 11 LIZ. f,t tit. yZAJxJO ^ 




V Z LKITI/SJ 


— DO.04U zt U.UUl 


zy.uz ( zt u.uui 


T[d] 


359.51 ± 0.09 


2.98565 ± 0.00003 


if [km/s] 


0.464 ± 0.003 


58.1 ± 0.4 


e 


0.847 ± 0.002 


0.013 ±0.013 


w[°] 


52.2 ±0.4 


11 ± 11 


to [JD-2.45X10 6 ] 


3998.1 ±0.1 


1497.5 ± 0.3 


mi [M©] 


1.43 


0.90 


m 2 sin / [Mjup] 


10.45 


0.38 


a [AU] 


0.9930 


0.03918 



parameters are T = 359.51 d, e = 0.847 and K = 464 m/s 
(Table 1). Details on the data analysis using radial velocities 
are given in Chapter 3. 

Assuming a stellar mass m\ = 1.43A'/0 (Tamuz et al, 
2008), Eqs. (25) and (66) can be used to derive a compan- 
ion minimum mass of sin / — 10.45 Mj up , orbiting the 
central star with a semi-major axis a = 0.99 AU. With the 
radial velocity technique it is impossible to determine the 
inclination /, and therefore we are unable to describe the 
orbit in three dimensions and to determine the exact mass 
of the planet. Nevertheless, the orbit in two dimensions 
(orbital plane) can be completely characterized. Astrom- 
etry is the only observational technique that can provide the 
full, three-dimensional orbit of the planet, but at present few 
planets have been observed by this method (Chapt. 6). 

In Figure 8b we have drawn the orbit of HD 156846 b. 
Because of its high eccentricity, the orbit is very elongated. 
As a consequence, the separation between the planet and the 
star ranges from 0. 15 AU at periapse to 1 .83 AU at apoapse. 
In our Solar System comets are the only objects that present 
such large variations in their position relative to the Sun. 
The origin of such high eccentricities is unknown, but a pos- 
sible explanation is through close encounters between very 
massive bodies during the formation process (Ford and Ra- 
sio, 2008). 

The angle to — 52.2° corresponds to the argument of 
periapse, that is measured from the nodal line between the 
plane of the sky and the orbital plane of the planet. For 
HD 156846 b this quantity is well defined, because the or- 
bit is so eccentric. According to Eq. (34), at periapse the 
orbital velocity is maximal and therefore it shows an easily 
identifiable peak in the observational data (Fig. 8). 

The planet is at periapse whenever 

t = t + kT with k = 0, ±1, ±2, ... , (67) 

where t = JD 2453998.1 (19 Sep. 2006 at 14h 24m UT) 
is the time of periapse passage. In fact, because the orbit is 
periodic, any instant of time, t , given by Eq. (67) can be 



used as the time of periapse passage. This is true for the 
two-body problem, but no longer valid if additional bodies 
are present in the system. Indeed, mutual planetary pertur- 
bations will disturb the orbits and the time of two succes- 
sive periapse passages is no longer given exactly given by 
Eq. (67) (see Chapt. 10). Although observers often use to 
as a parameter to characterize orbits, for multi-planet sys- 
tems it is meaningless. A better option is to use the mean 
anomaly M (Eq. 40) or the mean longitude A (Eq. 46). 
The observer fixes a date tf and then provides the value 
of M = M computed for that date (Eq. 40): 

M =n(tf-to). (68) 

A practical choice of tf is to use tf = t , since M = 
(and A = w). For multi-planetary systems, the planet 
will still be at periastron whenever M = 0, but the time 
of successive periapse passages will no longer be given by 
Eq.(67). 

6.2. HD 83443 b 

HD 83443 b is a short period Jupiter-size planet, there- 
fore belonging to the class of "Hot-Jupiters". It was first an- 
nounced at the Manchester IAU Symp. 202 as a resonant 2- 
planet system with periods T x = 2.986 d and T 2 = 29.85 d, 
but subsequent observations could not confirm the presence 
of the companion around 30 d. The origin of the transient 
signal is not clear yet, but an appealing possibility is to at- 
tribute the effect to activity of the star (Mayor et al, 2004). 
As a consequence, the HD 83443 star has been monitored 
many times. After 257 radial velocity measurements us- 
ing the CORALIE spectrograph with a mean uncertainty of 
8.9 m/s (taken from March 1999 to March 2003), the short- 
period planet at T = 2.986 d was found to revolve alone in 
an almost circular orbit (e <~ 0) with K = 58.1 km/s (Ta- 
ble 1). In Fig. 9a we show the CORALIE radial velocities 
and the corresponding best-fit Keplerian model. 

Assuming a stellar mass mi = 0.90M Q (Mayor et ah, 
2004), a companion minimum mass m 2 sin / = 0.38 Mj up 
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Fig. 8. — (a) Radial-velocity measurements as a function of Ju- 
lian Date obtained with CORALIE for HD 156846, superimposed 
on the best Keplerian planetary solution (Table 1). (b) Keplerian 
orbit of HD 156846 b and reference angles. 



and semi-major axis a = 0.039 AU has been derived. 
Again, it is impossible to determine the inclination I, and 
the orbit can only be characterized in terms of its orbital 
plane (Fig. 9b). 

Because the orbital eccentricity is small and uncertain 
(e = 0.013 ± 0.013), so too is the argument of the peri- 
apse (pj = 11° ± 11°). Indeed, for perfect circular orbits 
(e = 0) the argument of the periapse is not defined since 
the distance from the planet to the star is constant. There- 
fore, it is also meaningless to provide the time of periapse 
passage (t = JD 2451497.5 ± 0.3, i.e. 15 Nov. 1999 at Oh 
0m UT). The fact that the error bar is only 0.3 d, suggests 
(erroneously) that it is small. However, since the orbital pe- 
riod of the planet is only 2.986 d, an uncertainty of ±0.3 d 
is equivalent to a 20% uncertainty in to- 

For circular orbits the parameters u> and to are not de- 
fined (because e = 0) and the same is also true for the 
mean anomaly M (angle between the periastron and the 
planet). However, the orbit of the planet is still well deter- 
mined (Fig. 9a) and we should be able to provide accurate 
positions of the planet on its orbit. The correct parame- 
ter for that purpose is the sum M + lo or the mean lon- 
gitude X = M + m (Fig. 9b). Indeed, fixing the date at 
t f = JD 2453000.0 we obtain A = 91° ± 1°, which has a 
relative error of about 0.3%, much better than the 100% of 
uncertainty in w. 
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Fig. 9. — (a) Phase-folded radial-velocity measurements ob- 
tained with CORALIE for HD 83443, superimposed on the best 
Keplerian planetary solution (Table 1). (b) Keplerian orbit of 
HD 83443 b and reference angles. 
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